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ABSTRACT 

We present the X-ray broadband power spectral density function (PSD) of the X-ray-luminous 
Seyfert IC 4329a, constructed from light curves obtained via Rossi X-ray Timing Explorer monitoring 
and an XMM-Newton observation. Modeling the 3-10 keV PSD using a broken power-law PSD shape, 
a break in power-law slope is significantly detected at a temporal frequency of 2.5^'y x 10 -6 Hz, which 
corresponds to a PSD break time scale Tb of 4.6l 2 3 days. Using the relation between Tb, black hole 
mass Mbh, and bolometric luminosity as quantified by McHardy and coworkers, we infer a black hole 
mass estimate of Mbh = l-3lo 3 x 10 s M Q and an accretion rate relative to Eddington of 0.21^q'5o 
for this source. Our estimate of Mbh is consistent with other estimates, including that derived by 
the relation between Mbh and stellar velocity dispersion. We also present PSDs for the 10-20 and 
20-40 keV bands; they lack sufficient temporal frequency coverage to reveal a significant break, but 
are consistent with the same PSD shape and break frequency as in the 3-10 keV band. 
Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies — galaxies: individual (IC 
4329a) 



1. INTRODUCTION 

The rapidly variable, aperiodic X-ray continuum emis- 
sion from Seyfert AGN has long supported the notion 
that the X-ray emission originates from within a few tens 
of Schwarzschild radii of the central supermassive black 
hole. In the late 1980's, EXOSAT probed X-ray variabil- 
ity on time scales of a few days or less, yielding power 
spectral density functions (PSDs) covering temporal fre- 
quencies down to ~ 10 -5 Hz. These PSDs established 
the scale- invariant "red-noise" nature, i.e., larger vari- 
ability amplitudes toward larger time scales, character- 
ized by a PSD which increases toward lower temporal 
frequencies. However, no characteristic time scale, such 
as a "break" in the PSD power-law slope, was detected 
(e.g., Lawrence et al. 1987; McHardy & Czerny 1987; 
Green, McHardy & Lehto 1993). 

The Rossi X-ray Timing Explorer (RXTE) was 
launched in 1995; its unique attributes (large effective 
area, high throughput, fast slewing, and flexible schedul- 
ing) permitted multi-time scale monitoring campaigns 
that probed X-ray variability on time scales from hours 
to years (temporal frequencies from 10 4 — 10~~ 8 Hz). 
The resulting Seyfert broadband PSDs yielded evidence 
for breaks at temporal frequencies fh with PSD power- 
law slopes breaking from ~— 2 to ~— 1 above and below 
/b, respectively; the break frequencies corresponded to 
break time scales Tb of a few days or less (e.g., Edelson 
& Nandra 1999; Pounds et al. 2001; Uttley, McHardy 
& Papadakis 2002 (hereafter U02); Markowitz et al. 
2003b (hereafter M03); Papadakis, Reig & Nandra 2003; 
McHardy et al. 2004, 2005; Uttley & McHardy 2005). 

It was noted by M03 that a small sample of Seyfert 
PSDs was consistent with a relation between Tb and 
black hole mass Mbh estimated via reverberation map- 
ping, in the sense that relatively larger-mass black holes 
are associated with PSDs with larger values of Tb, albeit 
with some scatter. Using a larger sample, McHardy et al. 
(2006) demonstrated that the scatter could be explained 



by an additional dependence of Tb on TBoi/TEdd, the 
accretion rate relative to Eddington, in that for a given 
Mbh, sources accreting at relatively higher rates have 
PSDs shifted toward higher /b (smaller values of Tb). 
In the absence of an accurate estimate of Mbh by other 
means, X-ray timing can thus be used to estimate Mbh, 
if Tb and bolometric luminosity TboI are both known. 

In this paper, we present the broadband 3-10 keV PSD 
of IC 4329a, an X-ray bright and variable Seyfert located 
at a redshift of z = 0.01605 (Willmer et al. 1991). IC 
4329a had been monitored in the optical with the intent 
of using reverberation mapping to estimate Mbh in this 
object, but the constraints were quite poor (Winge et al. 
1996; Peterson et al. 2004). Here, we demonstrate that 
there is significant evidence for a break in the X-ray PSD 
of IC 4329a. We also present PSDs for the 3-5, 5-10, 10- 
20, and 20-40 keV bands. Section 2 reviews the multi- 
time scale sampling strategy and describes the observa- 
tions, data reduction, and light curve sampling. Section 
3 describes the methods of measurement and modeling of 
the PSDs. The results of model fits to the PSDs are pre- 
sented in §4. In §5, the PSD break frequency is used to 
yield an estimate of Mbh based on X-ray timing, which 
we compare to estimates obtained by various other meth- 
ods. 

2. OBSERVING STRATEGY, OBSERVATIONS, AND DATA 
REDUCTION 

We followed the observing strategy of Edelson and 
Nandra (1999), U02, M03, etc., in which a source is moni- 
tored with even, regular sampling on multiple time scales 
such that the resulting set of light curves yields individ- 
ual PSDs covering complementary temporal frequency 
ranges. With the reasonable assumption that the under- 
lying PSD shape has remained constant throughout all 
the monitoring observations 1 , the individual PSDs are 
combined to produce the final broadband PSD. 

1 The X-ray PSDs of black hole X-ray Binaries are observed to 
vary (display strong non-stationarity) on time scales of hours to 
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IC 4329a was monitored with RXTE once every 4.26 
days (64 orbits) for a duration of 4.3 years, from 2003 
Apr 8 to 2007 Aug 7 (Modified Julian Day (MJD) 
52737-54319; observation identifiers 80152-03-*, 90154- 
01-*, 91138-01-*, and 92108-01-*). This sampling, which 
probes variability on time scales from ~ a week to a few 
years, is henceforth called "long-term" sampling. Each 
observation lasted approximately 1 ks. There were four 
gaps due to sun-angle constraints in October-November 
of each year; each gap was ~ 50 days long. More inten- 
sive monitoring with RXTE was done to probe variability 
on time scales from several hours to a month ( "medium- 
term" sampling). RXTE observed IC 4329a once every 
~ 17.1 ks (three orbits) for a duration of 34.1 days, from 
2003 Jul 10 at 18:05 UT to 2003 Aug 13 at 20:35 UT 
(MJD 52830.75-52864.86; observation identifiers 80152- 
04-*). Again, each visit lasted ~ 1 ks. Finally, a contin- 
uous, 133 ks observation of IC 4329a with XMM-Newton 
(§2.3) allowed us to probe variability on time scales from 
~an hour to 1 day. 

2.1. PC A data reduction 

The Proportional Counter Array (PCA; Swank 
1998) consists of five large-area, collimated propor- 
tional counter units (PCUs). Reduction of the PCA 
data followed standard extraction and screening proce- 
dures, using HEASOFT version 6.5.1 software. PCA 
STAND ARD-2 data were collected from PCU 2 only: 
PCUs 1, 3 and 4 have been known to suffer from re- 
peated breakdown during on-source time, and PCU 
lost its propane veto layer in 2000. Data were rejected if 
they were gathered less than 10° from the Earth's limb, 
within 30 minutes of satellite passage through the South 
Atlantic Anomaly (SAA), if the satellite's pointing offset 
was greater than 0?2, or if ELECTRON2 was > 0.1). As 
the PCA is a non-imaging instrument, the background 
was estimated using the "L7-240" background models, 
appropriate for faint sources; see e.g., Markowitz, Edel- 
son & Vaughan (2003a) and Edelson & Nandra (1999) for 
details on PCA background subtraction, the dominant 
source of systematic uncertainty (e.g., in total broad- 
band count rate) in these data. Spectral fitting for each 
observation was done using XSPEC version 11.3.2ag, as- 
suming a Galactic column of 4.6 x 10 20 cm -2 (Kalberla 
et al. 2005). As the response of the PCA slowly hardens 
slightly over time due to the gradual leak of xenon gas 
into the propane layer in each PCU, 2 response files were 
generated for each separate observation. Fluxed light 
curves were extracted in the 3-10 keV band by fitting 
a power-law to that band only (2-10 keV is a "stan- 
dard" hard X-ray measurement band, but the 2-3 keV 
effective area of the PCA is extremely small, so we use 
3-10 keV here). Light curves for the 3-5 and 5-10 keV 
sub-bands and the 10-20 keV band were also extracted. 
Errors on each light curve point were derived from the 
standard error of 16-s count rate light curve bins within 
each observation. The total number of data points after 
screening was 372 points for each long-term light curve, 

days and longer. Scaling with Mbh and/or X-ray luminosity, we 
expect significant changes in the shapes of AGN PSDs only on 
timescales of centuries to millennia. 

2 In two previous papers, Markowitz, Reeves & Braito (2006) 
and Markowitz et al. (2009), the cause of the evolution in response 
was incorrectly listed as propane leaking into the xenon layers. 



with 14.5% of the data missing due to, e.g., sun-angle 
constraints or screening. Each medium-term light curve 
contained 173 points, with only 4.1% of the data missing. 
The 3-10 keV light curves are plotted in Figure 1. 

We also extracted 20-40 keV PCA light curves, as it 
would have been desirable to have PCA data simulta- 
neous with the 20-40 keV High-Energy X-Ray Timing 
Experiment (HEXTE) light curves (see below), but in 
this band, the source was only ~7% of the total PCA 
background; ±2% systematic uncertainties in the PCA 
background on time scales of ~l-2 ks ( Jahoda et al. 2006, 
Figure 29) yield uncertainties of ~ ±30% uncertainties 
in the net source count rate, and so we do not consider 
the PCA above 20 keV. 

2.2. HEXTE data reduction 

The High-Energy X-Ray Timing Experiment 
(HEXTE) aboard RXTE consists of two inde- 
pendent clusters (A and B), each containing four 
NaI(Tl)/CsI(Na) phoswich scintillation counters (see 
Rothschild et al. 1998) which share a common 1° FWHM 
field of view. Each of the eight detectors has a net open 
area of about 200 cm 2 . Source and background spectra 
were extracted from each individual RXTE visit using 
Science Event data and standard extraction procedures. 
The same good time intervals used for the PCA data 
(e.g., Earth elevation and SAA passage screening) were 
applied to the HEXTE data. To measure real-time 
background measurements, the two HEXTE clusters 
each undergo two-sided rocking to offset positions, in 
this case, to 1.5° off-source, switching every 32 s. There 
is a galaxy cluster (Abell 3571, located at a redshift of z 
= 0.039) located about 2° south of IC 4329a, at RA. = 
13h 47.5m, decl. = -32° 52 m. This source is detected 
in the RXTE all-sky slew survey (XSS; Revnivtsev et 
al. 2004), which shows the 8-20 keV flux of this source 
to be about half that of IC 4329a (at 2° off-axis, the 
count rate in the PCA is negligible, and thus the effect 
of contamination from Abell 3571 on PCA light curves 
and PSDs is negligible; Jahoda et al. 2006). However, 
BeppoSAX-PDS observations have shown no detection 
of Abell 3571 above 15 keV (Nevalainen et al. 2004). 
The 8-20 keV emission seen by the XSS must therefore 
be emission only between 8 and < 15 keV, and in this 
paper we use the 20-40 keV band, so the presence of 
A3571 can be safely ignored as far as contaminating 
HEXTE background data obtained within ~l-2° of the 
center of A3571 is concerned. Cluster A data taken 
during the following times were excluded, as the cluster 
did not rock on/off-source: 2004 Dec 13 - 2005 Jan 14, 

2005 Dec 12 - 2006 Jan 4, during 2006 Jan 25, and after 

2006 Mar 14. Detector 2 aboard cluster B lost spectral 
capabilities in 1996; these data were excluded. Light 
curves were extracted in 16 s bins in the 20-40 keV band. 
Deadtime corrections were applied to account for cluster 
rocking, pulse analyzer electronics, and the recovery 
time following scintillation pulses caused by high-energy 
charged particles; typically, the HEXTE deadtime is 
~ 30 — 40%. The light curves were then binned to every 
12.79 days (long-term) or 50.2 ks (medium-term). This 
action minimized variability associated with background 
systematics. The HEXTE background is relatively 
stable over long time scales but rapidly variable within 
a single satellite orbit as the spacecraft moves in the 
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geomagnetic environment. For a typical background 
rate, the Poisson error in a ^600 s good time exposure is 
>1% (Gruber et al. 1996). In the case of IC 4329a, the 
net source flux is roughly 9% of the total background in 
the 20-40 keV band, and each observation was only 1 ks 
in duration, yielding a >10% systematic uncertainty in 
the net source count rate. In addition, there are likely 
systematic uncertainties associated with the deadtimc 
correction, <0.5% in a 1 ks duration observation 3 . 
Errors on each point were determined by the standard 
error of the 16 s points. For observations using both 
clusters, light curves from clusters A and B were added; 
all light curves are in units of ct s _1 detector -1 (that 
is, ct s _1 per 7 detectors when both clusters were in 
operation, or per 3 detectors when only cluster B was 
operating normally; there were no observations where 
cluster A was the only cluster operating normally). The 
long-term light curve had 125 points, with 13 points 
(10.4%) missing; the medium-term light curve had 60 
points, with no points missing. 

The 20-40 keV light curves are plotted in Figure 1. 
The HEXTE gain and response are both very stable over 
time scales of years, permitting us to work in units of 
count rates. However, the HEXTE light curves, plotted 
in Figure 1, have had count rates converted to fluxes for 
consistency with the way the PCA light curves are plot- 
ted; we use a conversion rate of 1 ct s _1 per HEXTE 
detector in 20-40 keV corresponding to 2.6 xl0~ 10 erg 
cm~ 2 s _1 for a source with a photon index of 1.7. Vi- 
sually, there appears to be some mild variability in the 
20-40 keV light curves not present in the 3-10 keV light 
curve, likely from the aforementioned systematic un- 
certainties associated with background subtraction and 
deadtime correction. However, we note that the overall 
shape of the resulting long- and medium-term combined 
PSD (see §4) is roughly consistent with that obtained for 
the 10-20 keV band at the temporal frequencies probed, 
and we are confident that the systematic effects associ- 
ated with HEXTE background subtraction and deadtime 
correction do not significantly affect our conclusions for 
the 20-40 keV PSD. 

2.3. XMM-Newton data reduction 

XMM-Newton observed the nucleus of IC 4329a dur- 
ing 2003 Aug 6-7 (observation identifier 0147440101), 
for a duration of 136 ks. In this paper, we only con- 
sider data collected using the European Photon Imag- 
ing Camera (EPIC) pn camera. Light curves were ex- 
tracted using XMM Science Analysis Software version 
7.1.0, XSELECT version 2.4a, and the latest calibra- 
tion files. Source photons were extracted from a cir- 
cular region of radius 40"; backgrounds were extracted 
from circles of identical size, centered ~3' away from the 
core. Hot, flickering, or bad pixels were excluded. We 
inspected the 10-13 keV pn background light curve for 
flares, and excluded periods when the 10-13 keV back- 
ground rate exceeded 0.1 ct s _1 . We extracted a light 
curve in the 3-10 keV band binned to 2000 s; variability 
at shorter time scales was dominated by Poisson noise. 
The light curve had 67 points, with 3 points, or 4.5%, 
missing due to the above background screening. It is 
plotted in Figure 1. Light curves for the 3-5 and 5-10 

3 See http://mamacass.ucsd.edu/hexte/status/hexte_deadtime.html 



keV sub-bands were also extracted, binned to 5000 s (27 
points). 

3. PSD MEASUREMENT 

Light curve sampling parameters, including mean net 
source and background count rates, for all light curves 
are listed in Table 1. Also listed in Table 1 are variabil- 
ity amplitudes F var (see Vaughan et al. 2003 for definition 
of .Fvar and its uncertainty) for all light curves. Previous 
works have noted that, from 2 to ~20 keV, F var is gener- 
ally observed to decrease with increasing photon energy 
(e.g., Markowitz, Edelson & Vaughan 2003a; Miniutti et 
al. 2007); possible explanations include the dilution of 
the observed variability of the coronal power-law com- 
ponent by the Compton reflection hump or "pivoting" of 
the power-law component. The values of F var for the 3-5, 
5-10, and 10-20 keV bands are consistent with this trend 
(in the case of the 20-40 keV light curve, it is not im- 
mediately obvious how much variability is intrinsic and 
how much is due to the systematic effects discussed in 
§2). Table 2 lists PSD measurement parameters. The 
3-10 keV band has the best temporal frequency cov- 
erage; this band and the 3-5 and 5-10 keV sub-bands 
each contain long-, medium- and short-term sampling. 
As XMM-Newton lacks coverage above 12 keV, there is 
no short-term sampling to complement the 10-20 keV 
and 20-40 keV PSDs. Nonetheless, the 20-40 keV PSD 
presented here is, to our knowledge, the first X-ray PSD 
of a Seyfert above 20 keV. 

PSDs were measured for each light curve separately, as 
described in §3.1. However, as described in §3.2, these 
"observed" PSDs suffer from measurement distortion ef- 
fects (aliasing and red-noise leak) and the errors are 
not well-determined. A model-dependent Monte Carlo 
method was used to assign proper errors and to recover 
the intrinsic, underlying PSD shape. These steps are 
reviewed only briefly in the next two subsections; for 
further details, the reader is referred to U02 and M03. 

3.1. Initial PSD construction 

Initial PSD construction closely followed §3.1 of M03. 
Light curves were linearly interpolated across gaps, 
though such gaps were rare. Each light curve's mean 
was subtracted. Periodograms were constructed using a 
Discrete Fourier Transform (e.g., Oppenheim & Shafer 
1975; see also §3.1 of M03 for definitions of Fourier fre- 
quencies / and the periodogram) . 

Following Papadakis & Lawrence (1993) and Vaughan 
(2005), the periodogram was logarithmically binned ev- 
ery factor of 1.6 in / (0.20 in the logarithm) to pro- 
duce the observed PSD P(f); the two lowest temporal 
frequency bins were widened to accommodate three pe- 
riodogram points. The ^5 lowest-temporal frequency 
bins in each individual PSD typically contained <15 pe- 
riodogram points, precluding us from assigning normal 
errors. 

The individual long-, medium-, and short-term PSDs 
were combined to yield the final, broadband observed 
PSDs P(f), which are plotted in Figure 2, as well as in 
the top panels of Figures 3, 4, and 5 in / x Pf space. 
No renormalization of the individual PSDs was done. 
The constant level of power due to Poisson noise was 
not subtracted from these PSDs, but instead modeled in 
the Monte Carlo analysis below. 
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3.2. Monte Carlo Simulations 

Visual inspection of Figure 2 supports, at first glance, 
the notion that all five PSDs are at least roughly con- 
sistent in broadband shape. However, to account for 
PSD measurement effects (the reader is referred to U02 
and M03 for descriptions of aliasing and red-noise leak, 
and how the method accounts for them) and to assign 
proper errors, we use a version of the Monte Carlo tech- 
nique "PSRESP" introduced by U02. We summarize this 
technique as follows: one first assumes an underlying 
broadband model PSD shape, and, for each individual 
observed light curve sampling pattern, one uses the al- 
gorithm of Timmer & K6nig (1995) to generate Atrial 
light curves. The simulated light curves are each resam- 
pled to match the observed light curve, then PSDs are 
calculated. The model average PSD P s i m (/)> which ac- 
counts for the distortion effects, is calculated from the 
Atrial PSDs; errors for each PSD bin are assigned based 
on the RMS spread of the individual simulated PSDs at 
a given temporal frequency bin. The constant level of 
power due to Poisson noise Pp sn (see Table 2) is added 
to the model as Pp sn = 2(/z + B)//i 2 , where \x and B are 
the source and background count rates, respectively; for 
non-continuously sampled light curves, this must be mul- 
tiplied by the ratio of AT samp , the average sampling time, 
to the average exposure time per observation. For the 
HEXTE light curves, the estimates of Pp sn were calcu- 
lated taking into consideration the total number of good 
detectors (e.g., times when cluster A data were not used). 

The goodness of fit is determined. A statistic Xdist 
which compares the observed PSD P(f) to the model 
average PSD P s i m (/), using the errors from the model 
(as opposed to using poorly-determined errors from the 
observed PSD), is calculated. The best-fitting normaliza- 
tion of the P s im(/) is found by renormalizing P s i m (/) un- 
til the value of observed Xdist ^ s minimized. 10,000 com- 
binations of the long-, medium-, and short-term PSDs 
Psim(/) are randomly selected to model the Xdist distri- 
bution, comparing P s i m (/) to P s im(/)- The probability 
that the model PSD can be rejected, R, is defined as the 
percentile of the 10,000 simulated Xdist values exceeded 
by the value of the observed Xdist ■ 

A range of underlying PSD model shapes is thus tested 
in this way, and the model with the lowest rejection prob- 
ability (highest likelihood of acceptance L = 1 — R) can 
be identified. 

4. PSD FIT RESULTS 

We tested PSD model shapes consisting of simple un- 
broken power-law models and singly-broken power-law 
models; the quality of each PSD precludes more complex 
shapes such as single or multiple Lorcntzian components 
or quasi-periodic oscillations routinely modeled in the 
PSDs of X-ray Binary systems. 

The unbroken power-law model was of the form P(f) = 
Ao(///o)~ Q , where a is the power-law slope and the nor- 
malization Ao is the PSD amplitude at /o, arbitrary cho- 
sen to be 10 -6 Hz. Pp sn is added to each simulated PSD 
but is not explicitly listed here since Pp sn is different for 
each PSD segment. The model was tested by stepping 
through a from 0.0 to 3.2 in increments of 0.1, each time 
with At r iai=300 simulations done to calculate P s i m (/)- 



The best-fitting models are plotted in the second panels 
of Figures 3, 4 and 5. The best-fitting values of a, Aq, 
and likelihood of acceptance £ U nbr are listed in Table 3. 
The errors on a correspond to values la above the rejec- 
tion probability P un br for the best-fit value on a Gaus- 
sian probability distribution; for example, if the best-fit 
model had P un b r =95.45% (2. Oct on a Gaussian proba- 
bility distribution), errors correspond to P un b r =99.73% 
(3. Oct). The errors on Aq were determined from the RMS 
spread of the 10 4 randomly selected sets of simulated 
PSDs. 

The high rejection probabilities for the 3-10, 3-5, and 
5-10 keV PSDs (each > 90%) and the residuals plotted 
in Figures 3b and 4b suggest that a more complex PSD 
model shape may be appropriate for these three PSDs. 
The 10-20 keV and 20-40 keV PSDs cover less dynamic 
range in temporal frequency than their <10 keV counter- 
parts due to the lack of short-term sampling; the rejec- 
tion probabilities are consequently much lower, <90%, 
along with smaller residuals (Figure 5b). 

To test for the presence of breaks in the PSD, we then 
tested a singly-broken PSD model shape of the form 

pm _ f ^i(///b)- Ql °, /</b 
UJ I ^i(///b)- ai ", />/b 

where the normalization A\ is the PSD amplitude at the 
break frequency /b, and ol\ and ay are the low- and 
high-frequency power law slopes, respectively, with the 
constraint a\ < Qhi- (We also tested a more slowly- 
bending PSD model of the form P(/) = (Ai/- Ql °)/((1 + 
///b)^ Qhi-Q!l0 \ but given the PSD quality, there were de- 
generacies between /b, a\ a , and cthi such that reasonable 
constraints on /b could not be attained.) The range of 
slopes tested was 0.0-3.2 in increments of 0.1. Break 
frequencies were tested in the log from -7.4 to -4.9 in 
increments of 0.1, corresponding to / — > 1.26/ in the 
linear sc ale. 100 simulated PSDs were used to determine 
P s i m (/). The best- fit model parameters, along with like- 
lihoods of acceptance Lbrkn, are listed in Table 4; resid- 
uals are plotted in the bottom panels of Figures 3 and 4 
for the 3-10, 3-5, and 5-10 keV PSDs. Figure 6 shows 
contour plots of ay versus /b for these three PSDs at 
the respective best- fit values of 0t\ o . We use the ratio 
of the likelihoods of acceptance ibrkn/i-unbr between the 
broken and unbroken power law model fits to establish 
if a break is significant (we use £brkn/£unbr at least ~ 
10). For the 3-10 keV PSD, we conclude that a break 
is significantly detected at / b = i0"( 5 - 6 (+ - 5 -°- 3 )) Hz = 
2.51 + r 72 x 10 6 Hz, which corresponds to a time scale 

of 4.6^°3 1 days. Breaks are also significantly detected in 
each 3-5 and 5-10 keV sub-band PSD, with best values 
of /b also near 10 -6 Hz. Because of the limited dynamic 
range in temporal frequency for the 10-20 and 20-40 keV 
PSDs, the improvement in fit when adding a break to the 
model is not high, with £brkn/£ U nbr ~ 2 — 3. Parameters 
for the broken power-law model are thus listed in paren- 
theses in Table 4 for these two PSDs and data/model 
residuals are not plotted. 

We next explored the possibility that all five PSDs 
could be consistent with the same PSD shape, namely 
that of the 3-10 keV PSD, with a low-frequency power 
law slope equal to -1.0, a high frequency slope near -2, 
and a break frequency ~ 1 — 2 x 10~ 6 Hz. We re-tested 
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the 3-5 and 5-10 keV PSDs with a broken power-law 
model, keeping a\ a fixed at 1.0. Acceptable fits were ob- 
tained, with best-fit values for /b and «hi consistent with 
those measured for the 3-10 keV PSD; results are listed 
in Table 4. The 10-20 and 20-40 keV PSDs only cover 
up to roughly 10~ 5 Hz, and there is not much "leverage" 
above ~ l0~ 5 - 5 ~ 6 Hz to cither confirm or refute the pres- 
ence of a break at those temporal frequencies. However, 
the best-fit power law slopes in the unbroken power law 
fits are flat, 1.1-1.2 (compared to 1.5-1.7 for the < 10 
keV PSDs, which extend up to roughly 10~ 4 Hz, well 
above the best- fit 3-10 keV PSD break frequency). We 
conclude that the 10-20 and 20-40 keV PSDs are at least 
roughly consistent with the presence of breaks at or near 
that found in the <10 keV PSDs, and with the same a\ Q , 
and that the >10 keV PSDs probe temporal frequencies 
primarily below the break frequency. 

5. DISCUSSION AND CONCLUSIONS 

5.1. A New Black Hole Mass Estimate for IC 4329a 
from X-ray Timing 

We have used complementary RXTE and XMM- 
Newton monitoring of the X-ray bright Seyfcrt AGN IC 
4329a to measure the 3-10 keV PSD, and we find sig- 
nificant evidence for a break in the modeled power-law 
slope at a temporal frequency of 2.5l 1 "_ 7 x 10~ 6 Hz, which 

corresponds to a break time scale Tb of 4.6i2° 3 1 days. 
Best-fit power law slopes above and below the break are 
ahi = 2.3^04 and a\ = l.Oj^g, respectively. For a dis- 
cussion of candidate physical mechanisms to explain the 
turnover, the reader is referred to, e.g., Edelson & Nan- 
dra (1999) and Arevalo & Uttley (2006). 

We can derive a new estimate for Mbh from the PSD 
break, using the empirical relation between Tb, bolomet- 
ric luminosity TboI, and Mbh quantified by McHardy 
et al. (2006), log(T b (days)) = 2.1 log(M BH /10 6 M Q ) - 
0.981og(T Bo i/10 44 ergs^ 1 ) -2.32. The average absorbed 
3-10 keV flux from the long-term PCA monitoring is 
11.18 x 10~ n erg cm~ 2 s _1 . To find the unabsorbed 
flux, we used a spectral model in XSPEC consisting of 
a power law with photon index 1.74 (Markowitz, Reeves 
& Braito 2006), absorbed by the Galactic column, and 
five layers of absorption modeled by Steenbrugge et al. 
(2005): a column of cold gas intrinsic to the host galaxy 
with column density Ah = 1.7 x 10 21 cm~ 2 , and four 
ionized absorbers, with column densities 1.3, 0.32, 6, 
and 2 xlO 21 cm~ 2 and ionization parameters log£ 4 = 
-1.37, 0.56, 1.92, and 2.70 erg cm s _1 , respectively. (A 
fifth ionized absorber quantified by Markowitz, Reeves 
& Braito 2006, with best-fit modeled values of log£ = 
3.73 erg cm s^ 1 and Ah = 1.4 x 10 22 cm -2 , is expected 
to produce only a narrow Fe K absorption feature and 
is ignored.) To model each ionized absorber, an XSTAR 
component was used (Bautista & Kallman 2001). Using 
this model, a luminosity distance of 78.6 Mpc (Mould 
et al. 2000), and assuming a cosmology with H = 70 
km s _1 Mpc -1 and A = 0.73, the unabsorbed, rest- 
frame 2-10 keV luminosity L2-10 is thus estimated to 
be 1.1 x 10 44 erg s _1 . Using Marconi et al. (2004, their 

4 The ionization parameter £ = L- lon n^ 1 r^^ , where Li on is the 
isotropic 1-1000 Ryd ionizing continuum luminosity, n e is the elec- 
tron number density, and r is the distance from the central contin- 
uum source to the absorbing gas. 



Figure 3b), Lboi/-^2-io = 30, and iboi = 3.3 x 10 45 erg 
s _1 . Using the relation of McHardy et al. (2006), we de- 
rive M B h = l-3to.3 x 10 8 M Q (a 7% uncertainty in the 
distance to IC 4329a (Mould et al. 2000) translates into 
an additional ~ 7% uncertainty on Mbh)- The accretion 
rate relative to Eddington, Tboi/ L^dd, is thus estimated 
to be 0.2li°;°6. 

This estimate of Mbh is in agreement with several es- 
timates obtained by other methods. Reverberation map- 
ping has been known to yield highly accurate estimates of 
Mbh, but only if the data are of sufficiently high quality. 
In the case of IC 4329a, the estimate of Mbh is poorly 
constrained due to relatively low-quality optical spectra, 
and formally only an upper limit, 9.90t 1 {;g|x 10 6 M© (see 
Peterson et al. 2004 for details). However, many other 
methods of estimating Mbh in IC 4329a yield estimates 
closer to ~ 10 8 M Q . This includes other methods based 
on estimating the distance -Rblr from the central contin- 
uum source to the Broad Line Region, such as the empir- 
ical relation noted by Kaspi et al. (2000) between -Rblr 
and the optical continuum luminosity, "recalibrated" by 
Vestergaard & Peterson (2006; their equation 5). Using 
the value of H/3 FWHM (RMS) from Wandel, Peterson 
& Malkan (1999), 5960 ± 2070 km s" 1 , and a value for 
AT A (5100A) of 1.64 ± 0.21 x 10 41 erg s" 1 from Kaspi et 
al. (2000), the Vestergaard & Peterson (2006) equation 
yields M BH = l^lj;* x 10® M Q . Vestergaard & Peterson 
(2006; their equation 6) also provide a recalibrated for- 
mula for estimates of Mbh based on H/3 luminosity and 
optical continuum luminosity. Winge et al. (1996) report 
a mean H/3 flux of 3.331 x 10~ 13 erg cm~ 2 s _1 ; assum- 
ing a luminosity distance of 78.6 ±5.5 Mpc (Mould et al. 
2000) yields an H/3 luminosity of 2.45j£^ x 10 41 erg s" 1 . 

The resulting estimate of Mbh is thus 6.8^43 x 1O 7 M . 
A virial estimate based on the photoionization method 
(Wandel, Peterson & Malkan 1999) is 2.2 x 10 7 M Q . 
The black hole mass can also be estimated via the well- 
known relation between stellar velocity dispersion ct* and 
Mbh- We use equation 19 of Tremaine et al. (2002), 
M BH - lO 8 - 13±o - O6 M x (a*/(200 km g- 1 )) 4 - 02 ^- 32 . 
Oliva et al. (1999) report cr* = 218 ± 20 km s" 1 or 
231 ± 20 km s" 1 using 1.59/Ltm Si or 1.62 /im CO(6,3) 
features, respectively. Using the average of these two 
values yields M BH = 2.171^ x 10 8 M Q . Finally, Niko- 
lajuk, Papadakis & Czerny (2004) suggested a prescrip- 
tion to estimate Mbh based on short-term (~ 1 day du- 
ration) X-ray variability amplitude measurements, pro- 
vided that the variability amplitude measured probes 
temporal frequencies above the PSD break and that the 
PSD power law slope at high temporal frequencies is -2; 
this assumption is consistent with the current PSD re- 
sults. Nikolajuk, Czerny & Papadakis (2004) estimated 
M BH = 1-24 x 10 8 M Q based on short-term RXTE light 
curves, and Markowitz, Reeves & Braito (2006) esti- 
mated M BH = 8.64 ± 0.60 x 10 7 M Q using a 2.5-12 keV 
light curve from the 2003 XMM-Newton pn observation. 

Several of the X-ray timing studies discussed in §1 have 
noted the remarkable similarity between Seyferts and 
black hole X-ray Binaries in terms of their broadband 
PSD shapes; furthermore, the Xb-MBH - £Boi relation ex- 
trapolates down to stellar-mass black holes over 6-7 or- 
ders of magnitude in both Mbh and X-ray luminosity. 
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However, it is not immediately obvious from the shape of 
IC 4329a's PSD and from the derived value of LBoi/iEdd 
if IC 4329a is an analog to a high/soft or low/hard state 
black hole X-ray Binary system. There is not enough 
temporal frequency coverage below the break in the PSD 
of IC 4329a to determine if the power-law slope of ~-l 
persists for several decades of temporal frequency below 
the break, which would identify it as an analog of Cyg 
X-l in the high/soft state (e.g., Axelsson, Borgonovo & 
Larsson 2005); monitoring covering temporal frequencies 
down to and below 10 9 Hz would be needed to deter- 
mine the ultra-low temporal frequency PSD shape of IC 
4329a. 

5.2. The energy- dependent PSD properties of IC J^329a 

The PSDs presented here are roughly consistent with 
same underlying broadband PSD shape at all energies 
probed: a break near 2 x 10~ 6 Hz, and power-law slopes 
of ~ —1.0 and ~ —2.0 below and above the break fre- 
quency, respectively. The 10-20 and 20-40 kcV PSDs are 
probing temporal frequencies primarily below the break 
frequency but are consistent with the presence of a break 
in the PSD up to 20-40 keV. 

Arevalo & Uttley (2006) noted that an increase in 
/b with photon energy is expected in a model incor- 
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porating inwardly propagating variations in the local 
mass accretion rate modifying the central X-ray emis- 
sion (Lyubarskii 1997, Kotov et al. 2001), with rela- 
tively harder X-ray bands associated with more centrally 
concentrated emissivity profiles. However, no significant 
breaks were detected in the 10-20 and 20-40 keV PSDs. 
The current 3-5 and 5-10 keV PSDs lack the necessary 
high temporal frequency resolution to discern any de- 
pendence of /b on photon energy. A much more dense 
temporal frequency sampling covering at least the range 
~ 10~ 7 — <~ 10~ 4 - 5 Hz would be required to adequately 
test any energy dependence of /b in this object. 
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TABLE 1 
Light Curve Sampling Parameters 











Mean Source 


Mean Source 


Mean Background 




Bandpass 


Time 






Flux 


Count Rate 


Count Rate 


^var 


(kcV) 


scale 


Instrument 


ATsamp 


(erg cm~ 2 s _1 ) 


ct s- 1 


ct s- 1 


(%) 


q in 
o— 1U 


Long 


v>r* A 


4. ZD a 


11 Ifi v 1(1 — H 

11. lo X 1U 


1 1 c;'? 
11. CM 


o no 

z.yy 


1 1 oq _L n no 
1 /.Zo ± U.Uo 




Medium 


v>r* A 
rlA 


1 ( . 1 KS 


o 70 \y in — 11 

y . ( z x iu 


1(1 11 

1U. 1Z 


z.yo 


q An _l_ n i *7 

o.oy ± u.i / 




onort 


pn 


oUUU s 


7 ai \y in — 11 
/.Do X 1U 


A Qfi 

4.oU 


n nQ 
U.Uo 


Q *70 1 n 1 i? 
O. / Z ± U. ID 


3-5 


Long 


PCA 


4.26 d 


4.22 X lO" 11 


4.94 


1.05 


18.61 ±0.12 




Medium 


PCA 


17.1 ks 


3.60 x 10- 11 


4.28 


1.05 


9.95 ± 0.22 




Short 


pn 


5000 s 


2.83 x lO" 11 


2.77 


0.02 


4.01 ± 0.21 


5-10 


Long 


PCA 


4.26 d 


6.96 x 10- 11 


6.63 


1.94 


16.47 ±0.11 




Medium 


PCA 


17.1 ks 


6.12 x 10" 11 


5.84 


1.94 


8.34 ±0.23 




Short 


pn 


5000 s 


4.77 x 10- 11 


2.06 


0.01 


3.17 ±0.27 


10-20 


Long 


PCA 


4.26 d 


8.22 x 10" 11 


2.33 


2.15 


15.07 ±0.21 




Medium 


PCA 


17.1 ks 


7.33 x lO" 11 


2.07 


2.14 


7.53 ± 0.53 


20-40 


Long 


HEXTE 


12.79 d 


8.38 x 10- 11 


0.32 


3.54 


17.03 ± 1.19 




Medium 


HEXTE 


50.2 ks 


8.07 x lO" 11 


0.31 


3.44 


12.82 ± 1.31 


Note. — 


ATsamp is the sampling interval. PCA count rates arc for one PCU. 


For the HEXTE light curves, count rates are per detector, 



taking into account the number of good detectors per cluster and observations where cluster A data were not used. PCA fluxes were 
determined from fitting simple power-law models over the given energy range. 



TABLE 2 

PSD Measurement Parameters 



Bandpass 


Time- 


Temporal Frequency 


Pp S n 


(kcV) 


scale 


Range Spanned (Hz) 


(Hz" 1 ) 


3-10 


Long 


1.26 x 10~ 8 - 1.06 x 10~ 6 


80.1 




Medium 


5.83 x 10~ 7 - 2.30 x 10" 5 


4.27 




Short 


1.29 x 10" 5 - 1.81 x 10~ 4 


0.41 


3-5 


Long 


1.26 x 10~ 8 - 1.06 x 10~ 6 


181 




Medium 


5.83 x 10~ 7 - 2.30 x 10~ 5 


9.72 




Short 


1.28 x 10~ 5 - 6.30 x 10~ 5 


0.72 


5-10 


Long 


1.26 x 10~ s - 1.06 x 10~ 6 


144 




Medium 


5.83 x 10~ 7 - 2.30 x 10" 5 


7.62 




Short 


1.28 x 10~ 5 - 6.30 x 10" 5 


0.98 


10-20 


Long 


1.26 x 10~ 8 - 1.06 x 10~ 6 


608 




Medium 


5.83 x 10~ 7 - 2.30 x 10~ 5 


32.8 


20-40 


Long 


1.25 x 10~ s - 3.50 x 10~ 7 


5940 




Medium 


5.75 x 10~ 7 - 7.41 x 10~ 6 


180 



Note. — -Pp sn is the power due to Poisson noise (see text for details). 



TABLE 3 

Unbroken Power Law Model Fits to PSDs 



Bandpass 
(kcV) 


a 


A 
(Hz- 1 ) 


^unbr 

(%) 


3-10 

3-5 

5-10 

10-20 

20-40 


i 7 +0.4 
L -'-0.3 

1 5 +0 - 4 
L -°-0.2 

i r+0.5 
L -°-0.2 
1.2 ±0.2 
1.1 ±0.3 


2.5 ±0.5 x 10 3 
3.0±°,' 4 x 10 3 
2.3i'o j x 10 3 

2.1 ±0.2 x 10 3 
5.5±°i x 10 3 


2.2 
9.4 
4.0 
17.5 
22.3 


Note. — Results from fitting the PSDs with unbroken power law model, 
of acceptance for this model, defined as one minus the rejection probability. 


Aq is the amplitude at / = 


10~ 6 Hz. Lunbr is the likelihood 
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TABLE 4 

Broken Power Law Model Fits to PSDs 



Bandpass 




fb 












(keV) 


"lo 


(Hz) 




"hi 


(Hz- 


- 1 ) 


^brkn 


^brkn/^unbr 


3-10 


1 0+ ' 4 


2.51±f ™ x 10" 


-6 


9 O+0.8 
^•"-0.4 


1.8 ±0.2 x 10 3 


87.8 


39.9 


3-5 


fi+ - 5 
u ' D -0.6* 


2.51±5"ii x 10 ~ 


-7 


i 7 +0.2 
L -< -0.3 


4 >s+ - 5 
4 -°-0.4 


x 10 4 


97.1 


10.3 




1.0 (fixed) 


7 - 94 ie 2 36 6 x 10 


-7 


-i 7 +0.5 
x ' '-0.2 


7.5 ±0.7 x 10 3 


90.4 


9.6 


5-10 


i 1+0.3 
J — L -0.2 


s.ietiH x io- 


-6 


2 0+ - 6 
z - u -0.3 


9.5 ±0.9 x 10 2 


82.5 


20.6 




1.0 (fixed) 


l-OOlo 75 X 10 ~ 


-6 




4 fi+ - 5 

4 - D -0.4 


x 10 3 


79.2 


19.8 


10-20 


(i-H°oi) 


(i-26i;;* 63 x io- 


" 5 ) 


(1.70) 


(3.8t°;f 


x 10 2 ) 


(29.4) 


(1.7) 


20-40 




(3.16i 28 # x 10 


- 7 ) 




(4.4t„;° 


x 10 3 ) 


(75.9) 


(3.4) 



Note. — Results from fitting the PSDs with a broken power law model. A\ is the amplitude at the break frequency /j,; ot\ and a^i 
are the power law slopes below and above /t,, respectively. -Lbrkn is the likelihood of acceptance for this model, defined as one minus the 
rejection probability. An asterisk (*) denotes that the parameter uncertainty pegged at the limit of the parameter space tested. Best- 
fit parameters for the 10—20 and 20—40 keV PSDs are in parentheses, as the broken power-law model fits do not provide a statistically 
significant improvement over the unbroken power-law model fits. For the 10-20 keV PSD, the best-fit break frequency is located near the 
edge of the temporal frequency range sampled, so no reasonable constraints on a^i are obtained. 
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Fig. 1. — Long-, medium-, and short-term light curves of IC 4329a. The top (middle) panel shows the long- (medium-) term 3—10 keV 
PCA and 20—40 keV HEXTE flux light curves (the latter curves renormalized to the mean of the PCA light curves). The lower panel shows 
the 3—10 keV short-term EPIC pn light curve binned to 2000 s and in units of EPIC pn ct s . 
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Fig. 2.— The observed PSDs for the 3-10 (black), 3-5 (red), 5-10 (blue), 10-20 (green), and 20-40 keV (orange) PSDs. Filled circles, 
open squares, and filled diamonds denote PSDs derived from long-, medium-, and short-term monitoring, respectively. 
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I a) 3-10 keV obsd PSD + best-fit broken PL model 
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Fig. 3. — The top panel shows the observed 3-10 keV PSD, plotted in / X P(f) space to visually emphasize the turnover. Filled 
circles, open squares, and filled diamonds denote PSDs derived from long-, medium-, and short-term monitoring, respectively. The solid 
line denotes the best-fit singly-broken power-law model folded through the sampling window (i.e., containing the distortion effects of PSD 
measurement and power due to Poisson noise). The middle and bottom panels show the residuals to the best- fit unbroken and broken 
power law models, respectively. 
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a) 3-5 (black) + 5-10 (red) keV obsd 
PSDs + best-fit broken PL models 



£ o 

o ° 
Dl, 



a 1 « 
a. 

e 

5 T 



T3 O 
X 



73 O 




1 — I I 1 1 1 I 1 — I 1 1 III 1 — I 1 1 III 1 — I I 1 1 1 l i 

b) Unbroken PL 



ftt] 1 — I I I 1 1 I 1 — I I I III 1 — I I I III 1 — I I I I i n 

c) Broken PL 



10 ° 10 ' 10 ° 10 1 

Temporal Frequency (Hz) 

Fig. 4.— Same as Figure 3, but for the 3-5 (black) and 5-10 (red) keV PSDs. 
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.a) 10-20 (black) + 20-40 (red) keV obsd 
PSDs + best-fit unbroken PL models 




++H 1 1 I I I I 1 1 1 1 1 I I I I 1 1 1 1 1 I I I I 1 1 

_b) Unbroken PL '_ 

t * , ' '• * ^ 



10"° 10"' 10"° 

Temporal Frequency (Hz) 



10" 



Fig. 5. — The top panel shows the observed 20—40 keV PSD; symbols are the same as in Figure 3. The solid line denotes the best-fit 
"distorted" unbroken power-law model. The bottom panel shows the residuals to the best-fit unbroken power law model. 
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Fig. 6. — Contour plots of Ohi versus /j, for the best-fit singly broken power-law models for the 3-10 (left panel), 3-5 (right panel; black 
contours) and 5—10 (right panel; red contours) keV PSDs. Each contour represents a slice through the parameter space tested at the best-fit 
value of oi (listed in Table 4). Solid, dashed, dot-dashed, and dotted contours denote 68, 90, 95.4, and 99.0% rejection probabilities, 
respectively. 



